Perturbative matching of the staggered four-fermion operators for e'/e * 



Weonjong Lee 

Los Alamos National Laboratory, Mail Stop B285, Los Alamos, NM 87545, USA 

(February 1, 2008) 

Using staggered fermions, we calculate the perturbative corrections to the bilinear and 
four-fermion operators that are used in the numerical study of weak matrix elements for e /e. 
We present results for one-loop matching coefficients between continuum operators, calcu- 
lated in the Naive Dimensional Regularization (NDR) scheme, and gauge invariant staggered 
fermion operators. Especially, we concentrate on Feynman diagrams of the current-current 
insertion type. We also present results for the tadpole improved operators. These results, 
combined with existing results for penguin diagrams, provide the complete one-loop renor- 
malization of the staggered four-fermion operators. Therefore, using our results, it is possible 
^ . to match a lattice calculation of K'^ — if" mixing and K — > n% decays to the continuum 

NDR results with all corrections of 0{g^) included. 
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I. INTRODUCTION 



^ , The neutral kaon system has been much studied since it was observed that the strangeness states (K^, 

I TC^) in the quark model mix to produce short-lived and long-lived kaons {Ks, K^) in nature. The discovery 

. of Kl TTTT decays revealed that the weak interaction violates CP symmetry. In nature, CP asymmetry 

' happens in two ways: indirect and direct CP violations. The dominant effect parameterized by e (the indirect 

. CP violation) comes from the fact that Nature's neutral kaon mass eigenstates are not CP-symmetric. The 

phenomenon that the weak interaction itself allows Kl tttt decays directly is referred to as "direct CP 
violation" and it is parameterized by e'. 

The recent results of Re(e'/e) announced by the KTeV and NA48 collaborations supported the 
existence of the direct CP violation. Even though the Standard Model (SM) of strong and weak interactions 
provides a straight-forward qualitative understanding of CP violation in terms of a single phase (5), pertur- 
" bation theory does not allow a reliable quantitative calculation of the size of CP violation. In particular, in 

. K —^ TTTT decays, the direct CP violation is a result of destructive interference from various contributions. 

Theoretically, these come from strong interaction effects on hadronic matrix elements of the effective four- 
quark operators. Since the energy involved in these decays is at a scale of around 500 MeV, non-perturbative 
' tools such as lattice QCD must be used to calculate them. 

H I There have been a number of attempts to calculate e'/e on the lattice using clover/ Wilson fermions 

staggered fermions and domain-wall fermions |^-||. In this paper, we study the method using staggered 
fermions ||^. Recently there has been a numerical study with controlled statistical errors in calculating 
the AI = 1/2 amplitude, using staggered fermions However, this numerical study suffers from two 
major difficulties. The first is the puzzling crux in calculating (tttt | Oi \ K) directly in Euclidean space 
[ p^ . Therefore, one calculates hadronic matrix elements of (tt | | K) and (0 | | K) on the lattice, 
and reconstructs (tttt \ Oi \ K) out of these amplitudes, using chiral perturbation at leading order [ pl] |. 
Correspondingly, this method has uncertainties from the neglected higher orders in chiral perturbation such 
as final state interactions. The second difficulty is that the complete set of one- loop matching formulae for 
gauge invariant staggered fermion operators was not available. Only part of one-loop matching formulae for 
gauge invariant Bk and a complete set of one-loop matching relationship for gauge non-invariant Landau 
gauge operators were available in Ref. |l^,|l^,|l3| . Hence, even though the numerical simulation has reached a 
point to have statistical control over "eye" diagrams |^ , it was not quite possible to do the complete one-loop 
matching of the gauge invariant lattice results to the continuum hadronic matrix elements. 
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In this paper, we present the one-loop matching formulae between the continuum operators, defined in the 
NDR scheme, and the lattice gauge invariant operators of phenomenological interest for the current-current 
insertion type (corresponding to the "eight" contraction type in Ref. [||). These results, combined with 

existing results for penguin diagrams in Ref. , allow matching lattice calculation of K'^ — K mixing and 
K ^ TiTT decays (including A/ = 1/2 amplitudes) to the continuum NDR scheme with all corrections of 
0(5^) included. 

This paper is organized as follows. In Sec. Owe explain our notation for staggered fermions, bilinears 
and four-fermion operators on the lattice. Sec. p^Il| is devoted to the full description of renormalization and 
matching of staggered fermion composite operators. We start by presenting Feynman rules, and carry out the 
renormalization for bilinear and four-fermion operators on the lattice. We explain the tadpole improvement 
procedure for gauge invariant staggered fermion operators. We present the matching formula between lattice 
a nd con tinuum operators as our final results. Finite renormalization constants are summarized in Tables 
^ - KXII . Useful Fierz identities are described in Appendix^ and^. We close with some conclusions. 



II. NOTATIONS AND TERMINOLOGY 

In this section we will specify our notation for the action, fermion fields and composite operators on the 
lattice. We use the same notation for the staggered fermion action as in Ref. : 



3 = 0^^ ^ ^ 77^ (n) (x{n)U ^{n)x{n + /i) - x{n + ^)Ul{n)x{n)) + mx{n)x{n) 



2a 



(1) 



where n — (711,712,^3,^4) is the lattice coordinate and r]^{n) — (— 1)"^^ t-nn-i. 

In order to construct the four spin component Dirac field, we adopt the coordinate space method suggested 
in Ref. |l6| . In this method, we interpret 16 staggered fermion fields (x) of each hypercube in the coordinate 
space as 4 Dirac spin and 4 flavor components. The 4 flavor Dirac field Q{y) is defined as 

0-(y) = — ^ ^(7a)"^x(2/a) (2) 

where a is the Dirac spin index, i is the flavor index and Nf is the number of degenerate flavors {Nf = 4). 
Define also 

yA = 2y + A, with A e {0, 1}^ (3) 

and 

7A=7^72^^73^-^7^ • (4) 

In order to construct bilinear operators, we use a 7 matrix basis for both spin and flavor matrices. Using 
the Dirac fields Q, we can express a general form of quark bilinear operators: 

Q""\yKis'^ ® ^]^))Q^-'iy) = -Lx{yA)hF¥^)ABx{yB) (5) 



Here, a, (3 are Dirac spin indices and i, j are flavor indices. The matrix (75 ® (,f) is defined as 

(7s ® £.f)ab = -^Tr(7]^7s7B7],) . (6) 

This interpretation is called the coordinate- space method. 

We can translate this coordinate space method directly into the momentum space method as follows: 

67W^)ab - E -^('^^'''''i^^^^hD^i-l)''-'' . (7) 



CD 
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Note that the momentum-space matrices (75 (g) ^f) proposed in Ref. are unitarily related to the 
coordinate-space matrices (75 ^f)- 

The continuum hmit of the staggered fermion action on the lattice corresponds to QCD with four de- 
generate flavors {Nf — 4) |l5|. fn this limit, 75 represents the Dirac spin matrix and £,f belongs to SU(4) 
flavor symm etry gr oup. There are numerous choices to transcribe the lattice operator for a given continuum 
operator ]r3| , po| , pl| . We adopt the following conventional choice of bilinear operator transcription: 

Osriy) = ■^^x{yA){ls ^F)ABx{yB) U{yA,VB) ■ (8) 

AB 



Here, UiijA, J/s) is a product of gauge links that makes Osf gauge invariant p3| , p^ , pi 19|. The link matrices 
U{yA,yB) are constructed by averaging over all of the shortest paths between i/a and ys [ |T3| , p"8|jT9| ] , such 
that Osf is as symmetric as possible: 

U{yA,yB) = j^^U{yA,yA + ^pi)--- 

p 

C/(yA + Api+Ap2+Ap3 , z/b) , (9) 

where P is an element of the permutation group (1234) and 

- {B^ - A^)fi . (10) 

There are two kinds of four-fermion operators which have different color contraction structure. The general 
form of color two-trace operators is 

[xi'{yA){iF^)ABxi"{yB)][xt{yc){i^W^')cDXdiyD)] Uab{yA,yB) uUycyo) (n) 

and the general form of color one-trace operators is 

[xiHyA){is ® S,F)ABxtiyB)][xt{yc){is' ® S.F')cDXd{yD)] UadiyA^yo) Ucbiycys) (12) 

where /i, /2, /a and /4 label the continuum flavor (for example, u, d ,s, • • •). Note that we have introduced 
flavor in two ways: first, the 4 degenerate flavors hidden in the hypercubic index A of the staggered fermion 
field X'^iyA), and second, the continuum flavors labeled by the index /. 

For later convenience, let us introduce a compact notation. We can represent the bilinear and four-fermion 
operators in terms of S (scalar), V (vector) , T (tensor), A (axial) and P (pseudo-scalar). For example, the 
bilinears can be expressed as 



1 

A.B 



[v^ ^S] = ^^ Y,ix{yA)h^. ® i)ABx{yB)]K{yA,yB) (13) 
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K ^ ^] = 2^i'xiyA)il,.5 (5)ABXiyB)]U{yA,yB) (14) 

A,B 

Here, the continuum flavor indices are dropped out. For the four-fermion operators, we use the same notation 
as the bilinears but need to distinguish between color one trace and color two trace operators. This notation 
is exemplified by the following: 

[s X p][p X s]ii = ^ myA){mr5)ABx{yBm{yc)h^)cDx{yD)] 

/ A,B,C,D 

■U{yA,yB)U{yc,yD) (15) 

[A^ X P][v^, X s]i = ^ [x(2/a)(7p5 ® £.5)ABx{yB)][x{yc){it. ® ihoxiyo)] 



Nj 

f A.B,C,D 



■U{yA,yD)U{yc,yB) (16) 



Here note that the sub-indices /, // represent the color one trace and color two trace operators respectively. 

There are two completely independent methods to transcribe the lattice operators using Fierz transfor- 
mation: one spin trace formalism and two spin trace formalism fl^ . In this paper, we choose two spin trace 
formalism to construct the lattice operators and it is also adopted for our numerical study on e'/e. 
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III. RENORMALIZATION AND MATCHING 



In this section, we will review the renormalization of the continuum operators and explain the matching 
procedure between the continuum and the lattice operators at one loop level. We use the formulation of 

Rcf. rajiiiiijgMiii. 



A. Feynman rules 

We use the same Feynman rules as in Ref. |p^ . There was a typo in the Feynman rule for the gluon 
propagator in Eq. (78) of Rcf. Therefore, we provide the correct Feynman rule here. 

D"{k) = (1 - , (17) 

E/3 ^ sin {^akp) ^ sin {^akp)]^ 

where a is a gauge fixing parameter in the general covariant gauge and superscripts /, J represent the color 
indices in the adjoint representation of SU(3). 

All the Feynman diagrams are calculated at zero quark mass, which induces infrared divergences. We 
regulate the infrared divergences by adding a mass term to the gluon propagators both on the lattice and in 
the continuum. This dependence on the gluon mass cancels out, when we match the lattice and continuum 
renormalized operators. 



B. Bilinear Operators 

Let us consider a bilinear operator, Osf with spin S and flavor F. In the MS scheme, there are three 
distinct sets of commutation rules regarding 75: Naive Dimensional Regularization (NDR), Dimensional 
Reduction (DRED, DREZ), and t'Hooft-Veltman prescription (HV) |||,||,|||. Here, we choose NDR, the 
most commonly used scheme. At one loop level, we can write a general form of the renormalized continuum 
operators as 

iO^sT^'^) = E (^ss'Sff' + -f-,zrS>F'){0's";'^'^) ■ (18) 

S'F' ^ 

where (• • •) represent amputated Green's functions between external quark states. Here, the superscript (i) 
{i € {0, 1, 2, • • •}) represents the number of loops and 



^SF^S'F' = ^SS'^FF' 



Ts\og{^) + C^°^'+Rs . (19) 

K 



The fj, is the renormalization scale and k is the gluon mass, which regulates the infra-red. Fs is the anomalous 
dimension matrix: 

rs-^(as-l) (20) 

j^lt^l^lsl^l^L ■ (21) 

Note that the vector and axial currents have vanishing anomalous dimensions. C^""* is a finite constant 
term which depends on the regularization and renormalization scheme. For the NDR scheme 23 1 , 

as = (4,1,0,1,4) , (22) 
(jCont ^ 0, ^, 0, y ) for 75 = (/, 7^, 0-^^, 7m5, 75) respectively. (23) 

The Rs terms contain logarithmic dependence on the external quark momentum due to anomalous dimension. 
They are universal {i.e. independent of the regularization and renormalization scheme). 
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On the lattice, there are eight Feynman diagrams (Figure 2 in Ref. [Q) contributing to bilinear operators 



at one loop level. Note that only half of the self-energy diagrams (Figure 2 (e), (f), (g) and (h) in Ref. |13|) 



contribute since the other half are absorbed in the wave- function renormalization of the external quark fields. 



The explicit analytic results of these diagrams are given in Appendix A of Ref. 1 14|. The renormalized lattice 



bilinear operators at one loop level can be expressed in terms of tree level lattice operators as 

(4^) 



S'F' 

where 

= ^Sss'Sff'Ts logiaK,) + C^f^s'F' + Sss'Sff'Rs ■ (25) 

Here, Ts is the same as in the continuum (see Eq. (pO|)) since the anomalous dimension at one loop is 
independent of the regularization and renormalization schemes. Cgp*.*g,p, is the finite constant term unique 
to the lattice regularization scheme. Numerical values of Cgp.^gip, are summarized in Tables || and In 
Ref. [ p^lJl^ , it is proven that the C/4(l) symmetry of the staggered fermion action insures that for any spin- 
flavor structure the two bilinear operators Xs(js ® £F)xd and XsilSb ® £,F5)Xd are renormalized identically 
regardless of their distance. Hence, in Tables || and we provide the Cgp^^ip, values for only one of these 
two bilinear operators. Our values of the Cgp.g,p, constants are related to the dj (i = SF and j = S'F') in 
Tables 6 and 7 of Ref. as follows: 



1 



Sss'Sff' ( Cf + Ts log(7r) - ClP^,,p, 



CLatt 
op. c/ pi — — L/_F 



'SF;S'F' 

where 



Ctj - ((cts - 1)(2 log(7r) + 1) + is) 



(26) 
(27) 



Cf = — - - . (28) 

Here, note that Cy is calculated for the DREZ scheme and 

^Cont,DREZ ^ (-^^ g, 0, y ) for 75 = (/, cr^"^' 7^5, Ts) respectively. (29) 

The values of ts for DREZ are the same as in Ref. . 

*S = (^,0, i,0, i) for7s = (1,7^,(7^1., 7^5, 75) respectively. (30) 

The Csp%ipi values in this paper are in agreement with those of Ref. within the statistical uncertainty 
of finite integrals obtained using the Monte Carlo integration method. However, our results (and those of 
"22j) for the scalar bilinears [S x X] {X £ {S,V,T,A,P}) and the mixing of the bilinear [V^ x V^] 
: V^] are in disagreement with those in Ref. while the rest of our results agree. We have 
presented the one-loop renormalized operators in terms of tree level operators both in the continuum and 
on the lattice. Now we need to relate the lattice and continuum operators. The matching condition is based 
upon the observation that at tree level, the continuum operators are, by construction, the same as lattice 
operators in the limit of the zero lattice spacing (a ^ 0). In other words, at tree level as a ^ 0, the lattice 
operators are related to the continuum operators as follows: 

^Latm ^ QConm ^ ^^^^ (3^^ 

Here, 0{a) represents terms of order a or higher, which are supposed to vanish as a — > 0. At one loop (and 
higher) order, the lattice and the continuum operators will differ. A carefully constructed mixture of lattice 
operators is needed to reproduce the desired contimmm operator. From Eqs. (^8|), (^), ([si]), we can connect 
the lattice operators with the continuum operators at one loop level via 
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Cont{l) 
SF 



E 

S'F' 

E 

S'F' 



Sss'Sff' 



JSS'OFF' 



(47r 



f yCont 
7^ \^SF:,S'F' 



yLatt 
'^SF:,S'F' 



{An)' 



r(5ss"^i=^F'rsln(/ia)- 



(47r 



{6ss'6ff'C^°''' 



r^Latt 
^SF;S'F' 



^^Latt{l) 



(32) 



This is the matching formula between the lattice and continuum bilinear operators. 

Lepage, Mackenzie, and Parisi observed that the tadpole diagrams cause the large difference between the 
bare gauge coupling gf){a) and the renormalized coupling gjjgifJ' = [ p6|j27| . They proposed a mean field 
method for removing the dominant effect of tadpole diagrams. This method of t adp ole improvement has 
been applied to quark mass |^ and composite operator renormalization |22yi2| , [l^ , [l3| extensively. Here, let 
us explain the tadpole improvement procedure for the bilinear operators. The first step in this procedure is 
the rescaling of the staggered fermion fields and gauge link fields p^j2^. 



X 
X 



uoX 



uoX 
uo 



(33) 



A gauge invariant choice of the mean field scaling factor uq is [E7|J2 



-Re(TrUu) 



1/4 



(34) 



Under this rescaling, the bilinear operators transform as 



latt{l) 
SF 



o 



latt{l) 
SF 



1 



Nf 

. 1-AsF 



V'(2/a)(7S ® S,F)ABlp{yB)U, 



AB 



Nf 



X{yA){lS ® iF)ABX{yB)U. 



AB 



■ o 



latt{l) 
SF 



(35) 



Here, Asf = E J - 

bilinear operator. The point is that the tadpole improved operators O^"^^'' are better representations of the 



e I 

continuum operators O gp 



Cont{i) ^ Now we can rewrite Eqs. (| 
terms of the tadpole improved operators. 



corresponds to the distance between the quark and antiquark fields of the 

ators O^"^^-* are better representations of the 
™), the one loop relationship on the lattice, in 



S'F' 



{6sS'S 



FF' 



E i^SS'^FF' 



S'F' 



~2 

5 yLatt \ 

2 ^SF;S'F'; ■ "0 



(4^) 

~2 _ 

y r^Latt \unLatt(0)\ 
-(j:^'^SF-S'F')\^S'F' I 



Latt{0)^ 



(36) 



where 



ryhatt 

^SF;S'F' 

CLatt 
SF-S'F' 



Sss'^FF'^i^s^ogiaK 

(1 



c 



^SF-S'F' 



Latt 
SF;S'F' 

2 



SsS'^FF'Rs ■ 



A ^(47r)2A A 
12 



(37) 
(38) 



Here, note that the gauge coupling constant is also tadpole improved: g — > g . Regarding the tadpole 
improvement of gauge coupling constants, refer to Refs. |^,^ 28 1 for details. Correspondingly, we can 
rewrite Eq. (|32|), the relationship between the continuum and lattice operators, as follows: 
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^Cont(l) ^ 
S F / J 



S'F' 



E 

S'F> 



c r \ ^ ( rrCont ryLatt 

dsS'OFF' + J^;^ yZgp.g.p, - ZsF:,S'F' 

Sss'Sff' + (-4^-)2 ^ss'Sff'^s \n{fia) + 



^^Lattil) 
Ug,p, 



TP- / 

IJSS' OFF' L-5 — Ugp.g,p, 



\2 



{An 

This is the tadpole-improved matching formula for the bilinear operators. 

C. Four Fermion Operators 

Let us consider a general form of the staggered four-fermion operators on the lattice: 



O^r"/'^ . (39) 



[x(75 «> ^f)x][x(7S' «> Cf')x] 

There are 16"* — 65536 four-fermion operators. Hence, the mixing matrices are naively 65536 x 65536. 
However, we are interested in only a subset of these huge mixing matrices. In this paper, we focus on 
operators of phenomenological interest, especially relevant to the direct CP violation, e'/e. 

On the lattice, gauge non-invariant four-fermion operators such as Landau-gauge operators mix with gauge 
non-invariant lower dimension bilinear operators. This requires additional non-perturbative subtractions and 
it is significantly more difhcult to extract the divergent mixing coefficients in a completely non-perturbative 
way. For this reason, we are not interested in gauge non-invariant operators. Therefore, we study gauge 
invariant operators in this paper, since the gauge invariance protects operators from mixing with lower 
dimension gauge non-invariant operators 

Let us first classify the Feynman diagrams into two categories. One is the current-current insertion type 
and the other is the penguin insertion type p^ , p^ . In the case of the penguin diagrams, the perturbative 
corrections and the matching formula for the gauge invariant operators are given at one loop in Ref. [ p^ . 
However, for the current-current diagrams, the perturbative matching formulae at one loop level are not 
available for the complete set of the gauge invariant staggered fermion operators of phenomenological interest 
(e'/e). Therefore, in this paper, we calculate perturbatively the current-current diagrams of these operators 
and obtain the corresponding matching formulae at one loop. 

We begin by reviewing the renormalization of the four-fermion operators for the current-current insertion 
type in the continuum. The continuum operators of phenomenological interest can be expressed in terms of 
the practical basis introduced in Ref. pO[ . 

01 = (^i7pL-02)(V'37M^^4) 

02 = (V'i7m^'02)('(/'37m-RV'4) 

03 = -2(^iLV'2)(?3^V'4) (40) 

where L = (1 — 75) and R = (1 + 75). There arc two more operators included in the practical basis for 
completeness in Ref. pO| |. 

Oi - (^iLV^2)(?3i^4) 

O5 = -^(1 + l^) 5^(?i7M7.i^2)(?37.7MiV'4) (41) 

The odd parity parts of these operators are not important for the following discussion, but they are included 
for generality. As in the bilinear operators, we choose NDR as our continuum renormalization scheme. In 
this scheme, there are many ways to extend the definition of four-fermion operators to n = 4 — 2e dimensions, 
which corresponds to different choices of evanescent operators. Here, we adopt the convention of Ref. [ pO[ . 
In this convention, the perturbative corrections of four-fermion operators are derived from those of bilinear 
operators by imposing the Fierz symmetry and by loosening the charge conjugation. Using this convention, 
we can express a general form of the renormalized continuum four-fermion operators as 
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o.^[S;;;^] (42) 



2 

Here, note that = SijZf^""'^ for i {1, 2, 3}, and Zf °"* is a 2 x 2 matrix defined as 

Zf°^* = Zv+aCu + Zy+ACb — {zs+p + 6)C'c 

= 2s+pa + zv+ACb - [zv+A + 6)4 (44) 
The ©4 and ©5 operators mix with each other. 



2Cont ^ ^^2ZT - zs+p)da + (-1)0, - ^{zg+P - 2zt + 2)C, 

Zg""' = {2zT)Ca + zs+pCb - ^{zs+P + 2zT + 4)C, (45) 
Here, the color mixing matrices C are 

e-K"o'?). e,.l{^!,). ,40, 

The renormalization coefficients zy+A, zs+p and zt can be expressed in terms of the bilinear corrections 
given in Eq. (|I^) as follows: 

zv+A = zv + ZA (47) 

zs+p = zs + zp (48) 



Fx log(^) + C$°"* + i?x 

K 



(49) 



Here, X S {S, V, T, A, P} represents the spin structure of bihnear operators and Fx, C^°"* are given in 
Eqs. ( po[ ) and ( p3|) respectively. Note that Rx terms will cancel out when we match between the lattice 
and continuum operators. For the NDR and DREZ schemes, zy ~ and za = 0, and so zy+A = 0. This 
simplifies the Z,p°"* further: 

Z^°"' = -{zs+P + 6)dc 
■^f = zs+pCb — &Cc 

= ZS+pCa - 6Ce . (50) 



^2 
yCont 



As a summary, we present the results of the renormalization of the continuum four-fermion operators in the 
practical basis as follows: 

Here, note that Tij — SijTi and (7^^°"* — SijC^""*^ for i £ {1, 2, 3}. The anomalous dimension matrices are 



Fi — Tv+A 



F2 — Ty-A 



-2 -6 
-6 +2 

-16 

f6 -2 

r -r _ , -2 +6 

i3-i2[P^S] - I +16 
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54 



4 f -5 +3 

3 A -3 +13 

4 / +7 +3 
'3 A +6 -2 

4 / -2 +6 
'3 ■ I +3 +7 



These are universal and independent of the renormaUzation scheme. For the NDR scheme, the finite constant 
matrices are 

-^Cont /^Cont 



/-^Uont /^uont ^^T- ( +11/6 11/2 

Oi -^v+A - Y^^i - I 11/2 +11/6 



O2 — L-y-A 



-23/3 -3 
-1/2 +1/6 



r<Cont _ I^Cont _ ( +1/6 —1/2 \ 

^3 - ^2[p-s\ +23/3 ) 



r^Cont 
O44 



r^Cont 



2 / +2 

3 ^ -3 +11 

1 / +4 +4 

2 \ +5 +1 



^Gont _ 1 / +1 +5 
^54 - 2 I +4 +4 



= +^ ( \, ) (53) 



2 / +11 -3 
3^0 +2 

For further details of other continuum renormaUzation schemes, refer to Ref. [^0| and references in it. 

This completes our review of the renormaUzation of the the four-fermion operators of the current-current 
insertion type in the continuum. We now describe renormaUzation of staggered four-fermion operators of 
the same type on the lattice. For the purpose of matching, it is convenient to work with the practical basis 
as in the continuum. Hence, we choose the same practical basis on the lattice. Using the compact notation 
introduced in Section |l^ in the two spin trace formalism [p^, we can express the practical basis as follows: 

Of"" = [V; X P\{V^, X P] + [A^ X P\{A^, X P\ 
O^-" = [y^ X P][l/^ X P] - [A^ X P][A^ X P] 
Qhatt ^ _2(^[S' X P][Sx P]-[Px P][P X P] 
O^"** ^[Sx P][S X P] + [P X P][P X P] 

^-{[SxP][SxP] + [PxP][PxP]-J2 [T,. X P] [T,, X P]) (54) 



Latt 



Here, note that we choose the same flavor structure as the pseudo-Goldstone mode of the conserved U{\)a 
symmetry. 

The SU{A) flavor symmetry of the staggered fermion action becomes exactly conserved only in the limit 
of a = 0. On the finite lattice of a ^ 0, only the discrete SW/^^diag symmetry and U{l)v ® U{1)a symmetry 
are conserved exactly, which rotates the spin and flavor simultaneously In calculating the perturbative 
corrections to the operators of Eq. (p^, the operators with different flavors mix with one another, since 
the SU{4) flavor symmetry is not conserved at a ^ 0. Therefore, we can express a general form of the 
renormalized staggered four-fermion operators on the lattice as follows: 
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(O. 



/■r\Latt 
Latt(l)^ 



Latt 



/-nLatt 



Latt{Q)^ 



(47r)2 



Tij log(aK) 



CLatt I D 



{O 



Latt(0)\ 



(55) 
(56) 



where O 



Latt{0) 



may include operators with spin or flavor structure not included in the practical basis of 



Eq. (p^). The anomalous dimension matrices Fy are scheme- independent and the same as in Eq. (p^). The 
Rij contains the external momentum dependence and will cancel out when we perform the matching between 
the continuum and lattice operators. The scheme-dependent finite constant matrix terms (7,^"** are given 



in Tables 



Ill-XXn. The analytic formula to calculate these finite constant terms are provided in Ref. Il4| 



Regarding Eq. (C2) of Ref. Q , we have obtained a slightly different analytical result for the corresponding 



diagram as follows: 



g3(e) ^ 



(4^) 



- 2Zoooo(7s (Xi Cf)cz3(7S' <^£.f')c'D' 



1 



(1 - a)-Zoooo(7s ® S,f)cd{is' £.f')c'D' 



(1 



+ (l-a)l 

_r(_l)M.(S+_F) 



^M■(S'+F') 



{iMS ® iMF)cD{lMS' ® S,MF')c'D' 



(-1) 



(S+F),,+M-{S+F) 



+ (-1) 



(S' + F')^ + M-(S' + F') 



(-1) 



{S+F)^ + (S+F)^+M-{S+F) 



+ (-1) 



{S' + F')^ + {S'+F'),, + M-{S'+F') 



X [liivMS ® S,tiuMF)cD{jtivMS' ^ ^tivMF')c'D' 



(57) 



The difference is localized only to the gauge-dependent terms proportional to (1 — a), the coefficient of the 
covariant gauge fixing. Using Eq. (^), we confirmed the gauge invariance by checking that the summation 
of all the gauge-dependent terms (i.e. oc (1 — a)) of all the one-loop Feynman diagrams vanishes for each 
spin and flavor structure of the gauge invariant four-fermion operators. 

We have checked the results of finite constant terms (7^^°** in three independent ways. First, we checked 
the conserved U{1)a symmetry. For example, the following operators are connected under the U{1)a trans- 
formation |13|: 



[V^, X P] [V^ X P] ""^^ [Vf, X P] [A^ X S] , [A^ X S] [V^ x P] , [A^ x S] [A^ x S] 



[A,, X P][A^ X P] 

[P X P][P X P] ^-^^ [P X P][S X S] 
[S X P][S X P] ^-^^ [S X P][P X S] 



[A,, X P] [y^ X S] , [y^ X S] [A^ X P] , [V^ X S] [V^ X S] 
[SxS][PxP], [SxS][SxS] 
[PxS][SxP], [PxS][PxS] 



(58) 



Hence, we confirmed that all the finite renormalization constants of the operators in the right-hand side of 
Eq. (pq) can be reproduced from those in the left-hand side, using a simple U{\)a transformation. Second, 
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we checked the gauge invariance. In other words, we confirmed that the summation of the gauge dependent 
terms proportional to (1 — a) vanishes for individual gauge invariant operators in the practical basis. Third, 
we checked the Fierz identities of Eq. (A6). Any color one trace operator can be expressed as a linear 



combination of color two trace operators, using Fierz transformation, which holds valid under renormalization 
on the lattice (but not in the NDR scheme in the continuum). We explain the details of these lattice Fierz 
identities in Appendix ^ and ^ 



The values of C,^"*' presented in Tables IniHXl are in agreement with those extracted from Table II of Ref. 



The values of C/<"" presented in Tables XI-XXII 



are new. 



We have presented renormalized operators in the continuum and on the lattice. We now describe the 
matching formula. As in the case of bilincars, the tree level lattice operators are equal to the tree level 
continuum operators in the limit of a = 0. 



QLatt{0) _ QCont(0) _|_ q^^-^ 



(59) 



Using Eqs. (|51|), (|5g) and (|59|), we can obtain the matching relationship between continuum operators and 
lattice operators at one loop level. 



■Cont{l) _ 



E 



9 



(An) 



9^ ( 

-r,,iog(H + 7^( 



(471 



a 



Latt 



o 



Latt{l) 



(60) 



This is the one-loop matching relation between the lattice and the continuum four-fermion operators. 

We now turn to the tadpole improvement of the staggered four-fermion operators. First, let us consider a 
general color two trace operator: 



1 



^11 = [x{yA){lS ® S.F)ABUABX{yB)] [X{yc){is ® CF)cDl^CDX{yD)] 



(61) 



Under the rescaling transformation of Eq. (B3^ , this operator transforms as 



O 



II 



On = 



(62) 



In case of color one trace operators, the distance between quark and anti-quark fields of the bilinears does 
not define the length of the gauge link. Therefore it is convenient to express the color one trace operator as 
a linear combination of color two trace operators, using Fierz transformations. 



(63) 



where Hij is the Fierz transformation matrix. Under the rescaling transformation of Eq. (t 



^11 



(64) 



where and A^,^^, correspond to distances between quark and anti-quark fields of the two bilinears in 
Ojj respectively. Using perturbation, we can rewrite 0\ as 



0] = Y^H,,0\,- 



(47r) 



(47r)2 12 



9 — A-' — 



^11 



(47r)2 12 



(47r)- 



(65) 
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where 
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jk 



^SF 



LAq, 



S'F' 



(66) 



Here, note that the tadpole improvement for color two trace operators are diagonal (oc Sij), while the tadpole 
improvement for color one trace operators allows mixing with off-diagonal operators. As a summary, let us 
express the tadpole-improved four-fermion operators in terms of unimproved operators as follows: 



^Latt 



Latt 



/nLatt 



Latt 



(An) 



Latt 



(67) 



where 



T ■ — 

^13 — 



M,j 
N,j 

(An)' 



Nji = —Si-, 



A' 



SF 



A* 

l^S'F' 



(68) 



(69) 



The values of My and Nij are tabulated in Tables [Il-XXIl. Using Eq. (pTj), one can rewrite the one-loop 
renormalization relationship of Eq. ( |5^ ) on the lattice in terms of tadpole-improved operators. 



-^Latt{l) 
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(47r)2 



Tij log(aK) -I- C. 



Latt 



- Tij + Rij {O 



?Lati(0)v 



(70) 



Here, note that the gauge coupling is also tadpole-improved. Correspondingly, we need to rewrite the 
one-loop matching formula of Eq. (60) in terms of tadpole-improved operators as follows: 



■Cont{\) 



E 



(47r) 



,f,l0g(H + 7|^(5g°"* 



c. 



Latt 



-:^Latt{l) 



(71) 



This is our final result for the tadpole-improved matching relationship of four-fermion operators. 



IV. CONCLUSION 



We have studied the one-loop perturbative matching formula between the lattice and continuum bilinear 
and four-fermion operators that are used in the numerical calculation of e'/e. Our main contribution is that 
the results presented in this paper, combined with the existing results for the penguin diagrams, make it 
possible to match a lattice calculation of K ^ tttt decays to the continuum NDR results with all the 0{g'^) 
corrections included. 

We have studied all five four-fermion operators in the practical basis. The first three operators (Oi, O2, 
and ©3) are relevant to the numerical evaluation of e'/e. These operators mix with a large number of off- 
diagonal operators. However, we found that the coefficients of the off-diagonal mixing are small compared 
with those of the diagonal mixing. After tadpole-improvement, the diagonal mixing terms are also under 
control except for the {[P x P][P x P])n term, which receives a large correction at one loop. Various 
attempts to reduce this large perturbative correction of the ([P x P][P y. P])ii term by improving the action 
and operators are under development |]3l|| . 
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APPENDIX A: FIERZ TRANSFORMATION (I) 



First, we define the notation for Fierz transformation: 

Ms = I® I, Mv = XI ® 7/., Mt = Y^ cr^^ ® a^,^, 

Ma = X 7m5 ^ 7m5i Mp = 75 (g) 75 



(Al) 



where ® represents the direct product and (7^,^ = ^[7^,71^] with signature defined by {7/^,7^} = 2(5^i^. We 
also introduce the following compact notation for spin indices: 



{Mv)c,f}-a'l3' = ^{lp,)c,(i ® (7M)a'/3' 



(A2) 



The same notation also applies to the flavor structure S^p ® £,f by switching 7 to Now we can express the 
Fierz transformation of Mx for X G {S, V, T, A, P} as follows: 



M : 



/ Ms\ 
Mv ' 
Mt 
Ma 

\Mp j 



Maf3'-a'f3 = FMap-a'fS' , {Mi)ap>-a'l3 = [Fij] {M j) afi-a' f3' 



(A3) 



where 



1 1-1-11 

4-2 0-2 -4 

-6 -2 -6 

-4 -2 0-2 4 
1-1-11 1 



(A4) 



The spin and flavor matrices of the staggered four-fermion operators are Fierz-transformed separately 
using the same relationship of Eq. (A3). Let us define Tx as the spin part and Sx as the flavor part. For 
example, 

[rA(8)Sp] = [A^ X P][A^ X P] 
[Ts®Sp] = [S-kP][S-kP] 

Using this notation, we can express a general form of the Fierz transformation as: 

rp l/l/4;/3/2 „ _ p p rp l/l/2;/3/4 

1/1/ 
ill 



rp ^ i/i/4;/3/2 _ _ p p rp ^ 1/1/2 1/3/4 



(A5) 



where /, // represent color one trace form and color two trace form respectively. Note that the negative 
sign on the right-hand side of Eq. (|A5| ) comes from the anti-commuting characteristics of the fermion fields. 
Using these relationships of Eqs. ( |A3D and (|A5|), we obtain the Fierz transformation of the practical basis 
defined in Eq. (p|). 





+ rA)(5 


5Sp 


/l/4;/3/2 


[(rv 


-Ta)(i 


5Sp 


/l/4;/3/2 


2[(r5 


-Tp)(i 


$Sp 


/l/4;/3/2 


[(rs 


+ rp)(5 


$Sp 


/l/4;/3/2 



i[(ry + r^) ® (5s - 5y - St + + Sp)]{/^^^^^^ 



4 

-^[(r^-Tp) 



(^5 — — ^' 



+ -[(ry - Fa) ® (Ss - Sy - St + 



^ \l/l/2;/3/4 
^ \l/l/2 1/3/4 
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- lii^s + Tp - Tt) ® EpY/f'-'f'f' = +^[{Ts + Tp) ® (S5 - - St + + Sp)]{/=^^^^^ 

(A6) 

Here, note that these Fierz identities hold valid under a transformation of switching color-trace indices 
(J <-+ II) between the left-hand and right-hand sides. Similarly, we can also derive another Fierz identity 
used to calculate Bk in the one spin trace formalism ||l^ : 

[{Tv + Ta) ® (Sp + Es)]{'^'-'^'^' = \[{rv + Ta) X [Es ~Et + Sp)]{/^^^^^^ (A7) 



APPENDIX B: FIERZ TRANSFORMATION (II) 



As shown in Eq. ( |Aq ) , color one trace operators can be expressed as a linear combination of color two trace 
operators. These Fierz identities hold under renormalization on the lattice. The finite integrals for color one 
trace operators are completely different from those for color two trace operators. Therefore, it will provide 



a non-trivial check of our calculation if the results in this paper satisfies the Fierz identities of Eq. (A6) 



Here, we describe how to check the Fierz identities step by step. First, we calculate the complete set of 



one-loop perturbative corrections for both color one trace and color two trace operators in Eq. (A6). Usually, 
the operators in the right-hand side of Eq. ( [A6| ) mix with hundreds of operators at one-loop level. In order 
to compare these results with the operators in the left-hand side, one needs to Fierz-transform the mixing 
operators in the right-hand side into the flavor combination (/i/4;/3/2) as in the left-hand side. For this 
purpose of large scale Fierz transformation, there exists a handy formula [ p^ . Using the following group 
orthogonality relationship for the gamma matrices: 



j'^hA)al: 



^hA)lx^S^x^pp, (Bl) 



one can derive a handy tool for the large scale Fierz transformation: 

hs '»^F)ABhs' '^^f')a'B' = Y^'^ilslK ^1£.f)ab'{7s'1k 'S) ^Uf')a'b ■ (B2) 
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Using this tool, we confirmed that the lattice results in this paper satisfy the Fierz identities of Eq. (|A^). 
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OSF 


Ts 


^Latt 


lOl 


8 


+55.628(1) 




8 


+14.855(1) 


1(E) 


8 
8 


-10.578(2) 
-29.965(3) 


1(E) is 


8 


-47.804(2) 


7m ® 1 
7m •X'Cm 
7m ® 
7m ^ Cm.- 
7m ® C.-P 
7m ® Cm5 

/ M ^i^o 

7m ® ?5 


oooooooo 


0.000(0) 
+19.705(1) 
-13.387(2) 
+4.527(1) 
-29.654(4) 
-45.998(2) 
-13.411(3) 
-30.010(3) 


7mi' ® 1 
7m.- ® Cm 
7m.' ® Cp 
7m.' ® Cm.- 
7m.' ® Cmp 

7m.- ® CpcT 


-8/3 
-8/3 
-8/3 
-8/3 
-8/3 
-8/3 


-14.620(3) 
-0.425(1) 

-29.670(4) 
+7.731(1) 

-14.199(2) 

-45.396(2) 



TABLE I. Renormalization constants of bilinear operators at one loop. Osf represents bilinear operators. Fs is 
the anomalous dimension defined in Eq. (|20[). Cg^p^.gp is the finite constant term defined in Eq. (|2^) for the diagonal 
case: S' = S and F' = F. 



Osf 




'-^SF:S'F' 


7m ® 


7m ® Cm 


-4.055(0) 


7m ® ?m5 


7m ® 


+0.862(0) 


7m 8) C'P 


7m ® Cm.' 


+1.982(0) 


7m C-P 


7m ® Cmp 


-1.982(0) 


7m.' ® Cp 
7m.' ® ?P 


7m.' ® Cm 
7m.' ^ C.' 


+0.902(0) 
+0.902(0) 



TABLE IL Renormalization constants of bilinear operators at one loop. Osf represents bilinear operators. Os'f' 
corresponds to the mixing bilinear operators in Eq. (p^). C''^p.g,p, is the finite constant term defined in Eq. ( p^ ) for 
the off-diagonal case: S' ^ S and F' ^ F. 
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^Latt{0) 


color trace 








[SxV^][SxV^] 


I 





-14.395(3) 


+ CTI 


[S X V^.][S X V^,] 


II 





-3.992(1) 





[S X V^][S X K] 


I 





-1.110(0) 





[S X Af,][S X 4^] 


I 





-0.313(0) 





[S X A^][S X A^,] 


II 





-1.219(0) 





[S X A^][S X A,] 


I 





-0.119(0) 





[V^ X S][V^ X S] 


I 





+0.342(0) 





[V^ X S][V^, X 5] 


II 





-0.539(0) 







I 





+1.397(0) 





[V/.t X T]!i.t'][V^ X 


II 





-2.170(0) 





[v;, X r„,][i/p X T^p] 


I 





+0.565(0) 





[v;, X r„,][i/p X T^p] 


II 





-0.499(0) 





[v;, X r„,][v^ X r,.p] 


I 





-0.051(0) 





[Vp X T^,][V^, X r,.p] 


II 





+0.154(0) 





[Vp X r,p][Fp X Tp,] 


I 





-0.051(0) 





[i/p X r,p][i^p X Tp,] 


II 





+0.154(0) 





[Vp X r,p][Vp X r„p] 


I 





+1.205(0) 





[Vp X T,p][Vp X r„p] 


II 





-1.972(0) 





[Vp X T,p][l/p X T,^] 


I 





+0.502(0) 





[V^ X T't/p][V^ X T^i/T^] 


II 





-0.307(0) 





[Vp X P][Vp X P] 


I 


+2 


-24.167(2) 


+2CTI 


[Vp X P][Vp X P] 


II 


-6 


-7.028(3) 





[r^ti^ X V^][T^i/ X V^] 


I 





+18.270(2) 


—cti 


X V^][X'^i/ X V^] 


II 





-4.544(1) 





[T^ti^ X V^][Tyi[/ X V[/] 


I 





-1.952(0) 





[T^ti^ X V^][T^[/ X Vi/] 


II 





+2.527(0) 





X V^][Xj[ii/ X Vp] 


I 





-0.223(0) 





[Pp. X Vp][Tp, X Vp] 


II 





+0.669(0) 





[Tp, X Vp][Tp, X Vp] 


I 





-0.223(0) 





[Ppi^ X Vp] [Ppi/ X Vp] 


II 





+0.669(0) 





[Ppi^ X Vp] [Ppi/ X Vp] 


I 





-15.294(2) 


+CTI 


[Ppi^ X V^] [Ppi/ X Vp] 


II 





-4.387(1) 





[Ppiv X V^] [Ppi/ X V/^] 


I 





+1.110(0) 






TABLE III. Renormalization constants of a color one trace four-fermion operator at one loop, (Of"")/ defined m 
Eq. (^). The Vij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The Cj^"** matrix is defined in Eq. (^^. The Mij matrix represents the tadpole-improvement defined in Eq. (36). 
cti = (47r)^/12. All the greek indices are summed under the condition of n ^ u ^ p ^ rj. 
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^Latt{0) 


color trace 










I 





+0.018(0) 





l^'fj.i' X j4^] [T^ij/ X -A^i] 


II 





-1.061(0) 





[Tf_LLy X yl^][T'^(iy X A,y] 


I 





-0.055(0) 





[TfiL/ X yl^][T'^jy X Aij] 


II 





-0.191(0) 





[Tf_Li/ X Afj,^ [^M^^ -^p] 


I 





+0.032(0) 





[7^11/ X -A^] [Tfj,i/ X -Ap] 


II 





-0.096(0) 





[r^ti^ X Ap] [r^jy X A^] 


I 





+0.032(0) 





[r^tiy X Ap] [r^jy X A^] 


II 





-0.096(0) 





[7^11/ X Ap] [r^iy X Ap] 


I 





+0.742(0) 





[Tf_iu X Ap] [r^iy X Ap] 


II 





-1.219(0) 





[^yui^ X Ap] [r^iy X A^^] 


I 





+0.119(0) 





X S-Ji^p X S] 


I 





+0.342(0) 





[^p X S-Ji^p X S] 


II 





-0.539(0) 





[Afx X X"^!/] [A^ X T/iiv] 


I 





+1.397(0) 





[A^ X Tyii/] [A^ X Tfxi/^ 


II 





-2.170(0) 





[Ap. X T'p,i/][Ap. X Tp.p] 


I 





+0.399(0) 





[A^i X T'^[/] [A^i X Ti/p] 


I 





-0.083(0) 





[Ap, X r^t/] [Ap, X Ti/p\ 


II 





+0.249(0) 







I 





-0.083(0) 





[A^ X T'i/p][A^ X Tpiy] 


II 





+0.249(0) 





[Ap, X T'iyp][A^ X Ti/p] 


I 





+1.336(0) 





[Ap, X Tjyp][Ap, X Tiyp\ 


II 





-2.367(0) 





[Ap, X X'jyp][Ap, X TfyTj] 


I 





+0.399(0) 





[A^ X P][A^ X P] 


I 


+2 


-24.768(2) 


+2CT/ 


[A,, X P][A^ X P] 


II 


-6 


-5.225(1) 





[P X V^][P X V^] 


I 





+17.108(2) 


— CTI 


[P X V^][P X V,,] 


II 





-4.149(1) 





[Px i;.][Px K] 


I 





+0.805(0) 





[Px v;][Px K] 


II 





+0.915(0) 





[P X A^][P X A^] 


I 





+1.073(0) 





[P X A^][P X A^] 


II 





-1.061(0) 





[P X A^][P X A.] 


I 





+0.183(0) 





[P X A^][P X A.] 


II 





-0.191(0) 






TABLE IV. Renormalization constants of a color one trace four-fermion operator at one loop, (Of"")/ defined m 
Eq. (^). The Tij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The Cj^"** matrix is defined in Eq. (^^. The Mij matrix represents the tadpole-improvement defined in Eq. (36). 
cti = (47r)^/12. All the greek indices are summed under the condition of n ^ u ^ p ^ rj. 



18 



^Latt{0) 


color trace 




/~iL att 




[SxV^][SxV^] 


I 





"5.581(1) 





[S X V^.][S X V^,] 


II 





+1.860(0) 





[S X V^][S X K] 


I 





+1.114(0) 





[S X V^,][S X K] 


II 





-0.371(0) 





[S X A^][S X A^] 


I 





-1.311(0) 





[S X A^][S X A^] 


II 





+0.437(0) 





[S X A^][S X A,] 


I 





-0.236(0) 





[S X A^][S X A,] 


II 





+0.079(0) 





[Vf, X S][V^ X 5] 


I 





-1.107(0) 





[V,, X S][V^, X S] 


II 





+0.369(0) 







I 





-1.106(0) 







II 





+0.369(0) 





[v;, X r„,][v^ X T^p] 


I 





+0.150(0) 





[V^ X T^Jv; X T^p] 


II 





-0.050(0) 





[v;, X T^,][v^ X r,p] 


I 





+0.403(0) 





[Vp X rp,][v; X r,.p] 


II 





-0.134(0) 





[Vp X r,p][v; X Tp,] 


I 





+0.403(0) 





[Vp X T,p][Vp X Tp,] 


II 





-0.134(0) 





[Vp X r,p][Vp X r„p] 


I 





-2.776(0) 





[i/p X r,p][i^p X r,p] 


II 





+0.925(0) 





[V^ X 2~'i/p][V^ X Tl/?^] 


I 





+0.150(0) 





[V/.t X Tl-'p] [V/j, X 2~fy?y] 


II 





-0.050(0) 





[Vp X P][T/p X P] 


I 


-6 


-4.500(0) 





[Vp X P][Vp X P] 


II 


+2 


-58.523(7) 


+4CT/ 


[Ppi^ X V^][Ppi/ X V^] 


I 





-2.767(0) 





[Ppi^ X V^][Ppi/ X V^j 


II 





+0.922(0) 





[Ppi> X Vp][Tpi/ X V[/] 


I 





-0.416(0) 





[Ppi> X V^][T^pi/ X V[/] 


II 





+0.139(0) 





[Ppi> X V^][Tpi/ X Vp] 


I 





+0.765(0) 





[Ppi^ X Vp][7piy X Vp] 


II 





-0.255(0) 





[Ppi^ X Vp][Ppi/ X Vp] 


I 





+0.765(0) 





[Ppi^ X V^][X'p[/ X Vp] 


II 





-0.255(0) 





[Ppi^ X V^] [^pi/ X Vp] 


I 





-5.968(1) 





[Ppi> X Vp] [T'pt/ X Vp] 


II 





+1.989(0) 





[Ppi> X Vp] [T^jy X Vr]] 


I 





+1.946(0) 





[Ti_Liy X Vp] [T^jy X V7] 


II 





-0.649(0) 






TABLE V. Renormalization constants of a color two trace four-fermion operator at one loop, (Of "")// defined in 
Eq. (^). The Tij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The C/j"** matrix is defined in Eq. (|5^). The Nij matrix represents the tadpole-improvement defined in Eq. (39). 
cti = (47r)^/12. All the greek indices are summed under the condition oi ^ ^ u ^ p ^ rj. 



19 





color trace 




^Latt 


AT 




I 





-0.772(0) 





[T^ij^ X j4^] [T^ij/ X -A^] 


II 





+0.257(0) 





X yl^][T'^(iy X Ai^\ 


I 





-0.045(0) 





X yl^][T'^iy X -A;/] 


II 





+0.015(0) 





X Ap] [!2~|fj,iy X -Ap] 


I 





-1.705(0) 





[7pi/ X Ap] [T'p./y X Ap] 


II 





+0.568(0) 





[r^j/ X Ap] [T'p.jy X At^] 


I 





-0.529(0) 





[Tprf X A.p] [T^p,!/ X Atj] 


II 





+0.176(0) 





[A^ X S]\A^, X S] 


I 





-0.478(0) 





[A^ X S]\A^, X S] 


II 





+0.159(0) 





[A^ X T'^[/] [Ap, X T'^ii/] 


I 





-1.106(0) 





[A^ X T'^[/] [Ap, X T"^!/] 


II 





+0.369(0) 





[Ap, X T^i/\ [Ap, X T'p,p] 


I 





+0.150(0) 





[A.^ X 'T'^i/][A^ X "J/ip] 


II 





-0.050(0) 





[A^ X Tt/pjlA^ X Tup\ 


I 





-2.776(0) 





[Ap. X T'i/p][Ap, X Ti/p] 


II 





+0.925(0) 





[Ap, X T'i/p][Ap, X Turi\ 


I 





-0.656(0) 





[Ap, X T'i/p][Ap, X T,^'q\ 


II 





+0.219(0) 





[A^ X P][A^ X P] 


I 


-6 


-4.500(0) 





[A^ X P][A^ X P] 


II 


+2 


+1.501(2) 





fp X v„\\p X y„i 


I 





-3.153(1) 





[Px v;][Px v„] 


II 





+1.051(0) 





[Px i;,][Px K] 


I 





+0.416(0) 





[Px V;][Px K] 


II 





-0.139(0) 





[P X A^][P X A^] 


I 





-0.377(0) 





[P X A^][P X A^] 


II 





+0.126(0) 





[P X A^][P X A,] 


I 





+0.045(0) 





[Px A^][Px A,] 


II 





-0.015(0) 






TABLE VI. Renormalization constants of a color two trace four-fermion operator at one loop, (Of""")// defined in 
Eq. (^). The Vij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The Cj^"*' matrix is defined in Eq. ([s^). The Nij matrix represents the tadpole-improvement defined in Eq. (39). 
cti = (47r)^/12. All the greek indices are summed under the condition oi ^ ^ v ^ p ^ r). 



20 



^Latt{0) 


color trace 




^Latt 




[SxV^][SxV^] 


I 





+23.853(3) 


—cti 


[S X V^.][S X V^,] 


II 





+0.839(1) 





[S X V^][S X K] 


I 





-0.139(0) 





[S X V^,][S X K] 


II 





+0.416(0) 





[S X A^][S X A^] 


I 





-1.443(0) 





[S X A^][S X A^] 


II 





-0.842(0) 





[S X A^][S X A,] 


I 





+0.015(0) 





[S X A^][S X A,] 


II 





-0.045(0) 





[Vf, X S][V^ X 5] 


I 





+0.515(0) 





[V,, X S][V^, X S] 


II 





+0.254(0) 







I 





-3.977(0) 







II 





+1.064(0) 





[v;, X r„,][v^ X T^p] 


I 





-0.216(0) 





[V^ X T^Jv; X T^p] 


II 





+0.648(0) 





[v;, X T^,][v^ X r,p] 


I 





-0.051(0) 





[Vp X rp,][v; X r,.p] 


II 





+0.154(0) 





[Vp X r,p][v; X Tp,] 


I 





-0.051(0) 





[Vp X T,p][Vp X Tp,] 


II 





+0.154(0) 





[Vp X r,p][Vp X r„p] 


I 





-3.486(0) 





[i/p X r,p][i^p X r,p] 


II 





-0.409(0) 





[V^ X 2~'i/p][V^ X Tl/?^] 


I 





+0.053(0) 





[V/.t X Tl-'p] [V/j, X 2~fy?y] 


II 





-0.158(0) 





[Vp X P][T/p X P] 


I 


+16 


-13.668(2) 


+2CT/ 


[Vp X P][Vp X P] 


II 





-2.529(3) 





[Ppi^ X V^][Ppi/ X V^] 


I 





-24.725(2) 


+CTI 


[Ppi^ X V^][Ppi/ X V^j 


II 





+1.778(1) 





[Ppi> X Vp][Tpi/ X V[/] 


I 





+0.981(0) 





[Ppi> X V^][T^pi/ X V[/] 


II 





-2.943(0) 





[Ppi> X V^][Tpi/ X Vp] 


I 





-0.223(0) 





[Ppi^ X Vp][7piy X Vp] 


II 





+0.669(0) 





[Ppi^ X Vp][Ppi/ X Vp] 


I 





-0.223(0) 





[Ppi^ X V^][X'p[/ X Vp] 


II 





+0.669(0) 





[Ppi^ X V^] [^pi/ X Vp] 


I 





-23.594(2) 


+CTI 


[Ppi> X Vp] [T'pt/ X Vp] 


II 





-1.621(1) 





[Ppi> X Vp] [T^jy X Vr]] 


I 





-0.139(0) 





[Ti_Liy X Vp] [T^jy X V7] 


II 





+0.416(0) 






TABLE VII. Renormalization constants of a color one trace four-fermion operator at one loop, (Oa"")/ defined m 
Eq. (^). The Vij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The Cj^"** matrix is defined in Eq. (^^. The Mij matrix represents the tadpole-improvement defined in Eq. (36). 
Cti = (47r)^/12. All the greek indices are summed under the condition oi ^ ^ u ^ p ^ rj. 



21 



^Latt{0) 


color trace 










I 





+1.820(0) 





l^'fj.i' X j4^] [T^ij/ X -A^i] 


II 





-0.289(0) 





[Tf_LLy X yl^][T'^(iy X -A,y] 


I 





+0.049(0) 





[T^ti^ X X -A;/] 


II 





-0.146(0) 





X Ayj,] [2~/J.iy X -Ap] 


I 





+0.032(0) 





[7^11/ X -A^] [Tfj,i/ X Ap] 


II 





-0.096(0) 





[r^ti^ X Ap] [r^jy X A^] 


I 





+0.032(0) 





[r^ti^ X Ap] [r^jy X A^] 


II 





-0.096(0) 





[T^ti^ X Ap] X Ap] 


I 





+1.575(0) 





[T^ti^ X Ap] X Ap] 


II 





+0.446(0) 





[^yui/ X Ap] [r^iy X A^^] 


I 





+0.015(0) 





[T^^ii^ X Ap] [r^j/ X Ar^] 


II 





-0.045(0) 





[^^ X S-liA^ X S] 


I 





-0.515(0) 





[A^ X 5] [yip X S] 


II 





-0.254(0) 





X Tp,i/][Ap, X Tp,iy] 


I 





+3.977(0) 





[Ap. X Tp,i/][Ap, X Tp,iy] 


II 





-1.064(0) 





[A^i X 2~'p,[/] [Ap, X Tf_i_p\ 


I 





+0.050(0) 





[Ap, X T^p,!/] [Ap, X 2~'p,p] 


II 





-0.150(0) 





[Ayi X T/ilv] [Ayi X Ti/p^ 


I 





-0.083(0) 





[A^ X 7^/Ai/] [A^ X 7^i/p] 


II 





+0.249(0) 





[Ap, X Tiyp][A^ X 2~pjy] 


I 





-0.083(0) 





[Ap, X Tjyp] [Ap, X Tpjy] 


II 





+0.249(0) 





[Ap, X X'jyp][Ap, X Tl/p] 


I 





+3.354(0) 





[A^ X TL-p] [A.^ X r^/p] 


II 





+0.804(0) 





[A^ X 2~i/p] [A^ X ^1/77] 


I 





+0.050(0) 





[A^ X Tiyp][A^ X T^i/T^] 


II 





-0.150(0) 





[yip xP][Ap xP] 


I 


-16 


+14.269(2) 


-2CTI 


[^p xP][Ap xP] 


II 





+0.725(1) 





[P X Vp][P X V„] 


I 





+24.465(2) 


—cti 


[P X V^,][P X Vp] 


II 





-0.997(1) 





[Pxl^p][PxK] 


I 





-0.166(0) 





[P x l^p][P X K] 


II 





+0.499(0) 





[P X Ap][P X Ap] 


I 





-1.952(0) 





[Px Ap][Px Ap] 


II 





+0.685(0) 





[Px Ap][Px A.] 


I 





-0.079(0) 





[Px Ap][Px A,] 


II 





+0.236(0) 






TABLE VIII. Renormalization constants of a color one trace four-fermion operator at one loop, (Oa"")/ defined m 
Eq. (^). The Tij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The Cj^"** matrix is defined in Eq. (^^. The Altj matrix represents the tadpole-improvement defined in Eq. (36). 
cti = (47r)^/12. All the greek indices are summed under the condition oi ^ ^ u ^ p ^ rj. 



22 





color trace 




^Latt 




[SxV^][SxV^] 


I 





+5.581(1) 





[S X V^.][S X V^,] 


II 





-1.860(0) 





[S X V^][S X K] 


I 





-1.114(0) 





[S X V^,][S X K] 


II 





+0.371(0) 





[S X A^][S X A^] 


I 





-1.311(0) 





[S X A^][S X A^] 


II 





+0.437(0) 





[S X A^][S X A,] 


I 





-0.236(0) 





[S X A^][S X A,] 


II 





+0.079(0) 





[Vf, X S][V^ X 5] 


I 





+1.107(0) 





[V,, X S][V^, X S] 


II 





-0.369(0) 





[v;, X r„,][yp X T^,] 


I 





-1.106(0) 





[v;, X r„,][i^p X T^,] 


II 





+0.369(0) 





[v;, X r„,][v^ X T^p] 


I 





+0.150(0) 





[V^ X T^Jv; X T^p] 


II 





-0.050(0) 





[v;, X T^,][v^ X r,p] 


I 





+0.403(0) 





[Vp X rp,][v; X r,.p] 


II 





-0.134(0) 





[Vp X r,p][v; X Tp,] 


I 





+0.403(0) 





[Vp X T,p][Vp X Tp,] 


II 





-0.134(0) 





[Vp X r,p][Vp X r„p] 


I 





-2.776(0) 





[i/p X r,p][i^p X r,p] 


II 





+0.925(0) 





[l/p X r,p][Vi X T„,] 


I 





+0.150(0) 





[Vp X np][Vi X n,,] 


II 





-0.050(0) 





[Vp X P][Vp X P] 


I 


+6 


+4.500(0) 





[Vp X P][Vp X P] 


II 


-2 


-61.523(7) 


+4CT/ 


[Pp. X Vp][Tp, X Vp] 


I 





-2.767(0) 





[Pp. X Vp][Pp, X Vp] 


II 





+0.922(0) 





[Pp. X Vp][Pp,. X K] 


I 





-0.416(0) 





[Pp. X Vp][Pp, X K] 


II 





+0.139(0) 





P,„ X V,,]\T,,„ X Vn] 


I 





+0.765(0) 





[Pp. X Vp][Pp. X Vp] 


II 





-0.255(0) 





[Pp. X Vp][Pp. X Vp] 


I 





+0.765(0) 





[Pp. X Vp][Pp. X Vp] 


II 





-0.255(0) 





[Pp. X Vp][Pp. X Vp] 


I 





-5.968(1) 





[Pp. X Vp][Pp. X Vp] 


II 





+1.989(0) 





[Pp. X Vp] [Pp[/ X Vy] 


I 





+1.946(0) 





[Pp. X Vp] [Pp[/ X V7] 


II 





-0.649(0) 






TABLE IX. Renormalization constants of a color two trace four-fermion operator at one loop, (O2"")// defined in 
Eq. (0). The Vij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (Bih. 
The Cj^"*' matrix is defined in Eq. (|5^). The Nij matrix represents the tadpole-improvement defined in Eq. (39). 
cti = (47r)^/12. All the greek indices are summed under the condition oi ^ ^ u ^ p ^ rj. 



23 



C9, 


color trace 




^Latt 




[Tfxi/ X [T^/j/ X A.f_i] 


I 





+0.772(0) 





[Tfxu X j4^] [T^ij/ X A.f_i] 


II 





-0.257(0) 





[Tf_LLy X yl^][T'^(iy X Aiy] 


I 





+0.045(0) 





[Tf_iu X yl^][T'^iy X -A;/] 


II 





-0.015(0) 





[Tf_Li/ X yip] [!2~|fj,iy X Ap\ 


I 





+1.705(0) 





[Tfiu X A.p\ [Tfxij X A.p] 


II 





-0.568(0) 





[Tfxi/ X yip] [r^jy X Ajj] 


I 





+0.529(0) 





[TjLiy X yip] [r^j/ X Arjl 


II 





-0.176(0) 





[A^ X S][A^ X S] 


I 





-0.478(0) 





X S][A^, X S] 


II 





+0.159(0) 





[yl^ X T'^t/] [yip, X Tf^ny] 


I 





+1.106(0) 





[yl^ X T'^[/] [yip, X T^^,] 


II 





-0.369(0) 





[yip, X T^i/] [yip, X T'p,p] 


I 





-0.150(0) 





[yl/i X 'T'^i/][yl^ X Tfxp^ 


II 





+0.050(0) 





X '71/p][yl^ X Ti/p\ 


I 





+2.776(0) 





[yip. X T'i/p][ylp, X Ti/p] 


II 





-0.925(0) 





[yip, X T'i/p][ylp, X T,yri] 


I 





+0.656(0) 





[yip, X T'i/p][ylp, X T,yri] 


II 





-0.219(0) 





[A^ X P][A^ X P] 


I 


-6 


-4.500(0) 





[^^ xP][A^ xP] 


II 


+2 


+1.498(2) 





\p X v„ifp X y„i 


I 





+3.153(1) 





[Px v;][Px v„] 


II 





-1.051(0) 





[Px i;,][Px K] 


I 





-0.416(0) 





[Px V;][Px K] 


II 





+0.139(0) 





[P X A^][P X A^] 


I 





-0.377(0) 





[P X A^][P X A^] 


II 





+0.126(0) 





[P X A^][P X A,] 


I 





+0.045(0) 





[PxA^][PxA,] 


II 





-0.015(0) 






TABLE X. Renormalization constants of a color two trace four-fermion operator at one loop, (©2 "")// defined in 
Eq. (j54|). The Vij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The C/j"** matrix is defined in Eq. ([s^). The Nij matrix represents the tadpole-improvement defined in Eq. (39). 
CTi = {A-k)'^ /\2. All the greek indices are summed under the condition oi ^ ^ v ^ p ^ r). 



24 



o 



color trace 



Mi- 



% X V^][V^ X V^] 
Vm X V^,][V^ X V^,] 
X, X X K] 

X, X X K] 

X, X K]X^ X Vp] 
X X X Vp] 



S X S-JIS X 5] 
5 X S] IS X 5] 

SxT^,.][SxT„,] 
S x T^„][S X Tf,,] 
S X T^,.][S X T,.p] 
S X T^„][S X Tpp] 
S X P][5' X P] 
S X pji^* X P] 



I 

II 
I 

II 

I 

II 

I 

II 










2 X (+2) 
2 X (-6) 



-0.155(0 
-0.507(0 
-1.068(0 
-0.458(0 
-0.349(0 
-0.150(0 
2 X (+25.186(2) 
2 X (+3.972(1) 










2 X (-2cTi) 




I 

II 

I 

II 

I 

II 



+33.748(1 
-2.258(1 
-0.581(1 
+0.194(0 
+1.249(0 
-0.416(0 











Vp X 
Vp X 
Vp X 
Vp X 
Vp X 
Vp X 
Vp X 
Vp X 

v;, X 

K X 



^m][^m X 
^mIXm X 
A^][Vp X 

^.]Xm X 
Ap][Vp X 

^.]Xm X 
A,][Vp X 

^.IXp X 

A,][Vp X 



A.] 

A,A 
A,] 
A.] 
Ap] 
Ap] 
Ap] 
Ao] 



I 

II 

I 

II 

I 

II 

I 

II 

I 

II 



-0.402(0 
+0.776(0 
+0.461(0 
+0.198(0 
-0.064(0 
+0.191(0 
-0.064(0 
+0.191(0 
+0.104(0 
+0.045(0 



Tjj,i/ X Ti_i 

Tjjii, X Tjj^ 

Tj_ii/ X Tjj^ 

Tj_ii/ X 

Tj_ii/ X Tj_i 

Tjjiu X Tjj^ 

Tjjiu X Tjj^ 

Tf_ii/ X Tj_i 

Tj^ii/ X 

Tf^ii/ X 

Tj^ii/ X Tp' 

71(1/ X 



[T^i/ X Tp' 



X 

X Tp 
X 
X 

X n 
X r„ 

X Tp. 
X Tp^ 
xTm, 
xTm, 

X To; 



I 

II 

I 

II 

I 

II 

I 

II 

I 

II 

I 

II 



+0.815(0 
-2.065(0 
+0.449(0 
-0.150(0 
-0.449(0 
+0.150(0 
+0.102(0 
-0.307(0 
+0.102(0 
-0.307(0 
-0.552(0 
+1.275(0 



TABLE XL Renormalization constants of a color one trace four-fermion operator at one loop, (Ol""**)/ defined in 
Eq. (js^l). The Tij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The C/j"** matrix is defined in Eq. (^^. The Mij matrix represents the tadpole-improvement defined in Eq. (36). 
cti = (47r)^/12. All the greek indices are summed under the condition oi ^ ^ v ^ p ^ i). 
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o 



Latt{0) 



color trace 



[A^ X 
[A^ X 

[^M X 

[A^ X 
[A^, X 
[A^ X 

[^M X 

[^M X 
[^M X 
[^M X 



V^][A^ 
V,][A, 

V.][A, 
V^][A^ 
V^][A^ 
V^][A, 

V^][A^ 



xV,] 
X K] 
X K] 
X K] 

X 

X V^] 
X V^] 

X v;,] 

X V-p] 



[Af, X A^] [^p X Af,] 
[Af, X A^i [^p X Af,] 
[Af, X A„] [A^ X A^] 
[A^ X A^][A^ X A^j 
[A^ X A^][Af, X Ap] 
[Ap X A^] [A^ X Ap] 



I 

II 
I 

II 
I 

II 
I 

II 
I 

II 



+32.138(1 
+2.572(1 
-0.452(1 
-0.194(0 
+0.305(0 
-0.915(0 
+0.305(0 
-0.915(0 
+0.971(0 
+0.416(0 



-2CTI 













I 

II 
I 

II 
I 

II 



+0.220(0 
-1.090(0 
+0.593(0 
-0.198(0 
+0.134(0 
-0.045(0 



[P X S][P X S] 
[P X S] [P X S] 



[P X T^,][P X T^p] 
[P X Tf,,][P X T^p] 
[P X P][P X P] 
[P X P] [P X P] 



I 

II 
I 

II 
I 

II 
I 

II 










x(-2) 
X (+6) 



2 X (- 



+0.155(0 
+0.507(0 
+1.068(0 
+0.458(0 
+0.349(0 
+0.150(0 
-39.872(2) 



2 X (+40.087(1) 










2 X (+2cTi) 




TABLE XII. Renormalization constants of a color one trace four fermion operator at one loop, (el"'")/ defined m 
Eq. (js^l). The Fij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The C/j"*' matrix is defined in Eq. (^^. The Mij matrix represents the tadpole-improvement defined in Eq. (36). 
Cti = (47r)^/12. All the greek indices are summed under the condition oi ^ u ^ p ^ rj. 
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o 



Latt{0) 



color trace 



[SxT^ 
[S X P][S 



S X T^,] 
S X Tp,] 
5 X T,.p] 
5 X Tf,.p] 
X P] 



I 

II 
I 

II 
II 








2 X (-16) 



-0.915(0 
+0.305(0 
-0.299(0 
+0.100(0 
2 X (+95.607(4) 



2 X ( 








-6CTl) 



[V^ X A,, 
[V^ X Afj, 
[V^ X 
[V^ X 
[V,, X 4p 
[V,! X 
[V,, X A^ 
[V^, X A^ 
[V^ X 4^ 
K X A^ 



Vf, X A^, 

Vj, X 

Vu X A^ 

Vf, X 
X Au 
X A„ 
X Ap 

Vp X Ap 

v;, X A, 
K X A. 



I 

II 
I 

II 
I 

II 
I 

II 
I 

II 



+0.682(0 
-0.227(0 
+0.395(0 
-0.132(0 
+0.191(0 
-0.064(0 
+0.191(0 
-0.064(0 
+0.089(0 
-0.030(0 



[A^ X V^, 
Up X 

[Ap X K 

[A^ X K 
Up X Vp 

[Ap X Vp 
[Ap X K 
[Ap X K 

[Ap X K 
[Ap X K 



Ap X Vp 
Ap X Vp 
Ap X K 
Ap X 14 
Ap X K 

Ap X K 
Ap X "l/p 

Ap X Vp 

Ap X "1/p 

Ap X Vo 



I 

II 
I 

II 
I 

II 
I 

II 
I 

II 



+6.015(1 
-2.005(0 
-0.387(1 
+0.129(0 
-1.530(1 
+0.510(0 
-1.530(1 
+0.510(0 
+0.833(0 
-0.278(0 



[PxTp 
[P X Tp 
[PxTp 
\P X T„ 



P X Tp 
P X Tp 
P X Tp 
P X 



[P X P][P X P] 



I 

II 
I 

II 
II 








2 X (+16) 



2 X ( 



+0.915(0 
-0.305(0 
+0.299(0 
-0.100(0 
+111.255(2) 



2 X (- 








'2CTl) 



TABLE XIII. Renormalization constants of a color two trace four-fermion operator at one loop, (CI3"")// defined 
in Eq. (^J). The Vij matrix represents an anomalous dimension defined in Eq. ( p^ and its values are given in 
Eq. " — 
Eq. 



The Cij" matrix is defined in Eq. (jS^)- The Nij matrix represents the tadpole- improvement defined in 
|69). cti = (47r)^/12. All the greek indices are summed under the condition oi ^ ^ v ^ p ^ i). 



27 





color trace 




/~i Lcitt 




[S X S][S X 5] 


I 





-0.077(0) 





[5 X S\[S X 5] 


II 





-0.254(0) 





[SxT^,.][SxT„,] 


I 





-0.977(0) 





[S x T^„][S X Tf,,] 


II 





-0.419(0) 





[S X T^,.][S X T,.p] 


I 





-0.099(0) 





[S X T^„][S X Tf,.p] 


II 





-0.042(0) 





[S X P] [5" X P] 


I 


-2 


-25.186(2) 


+2CT/ 


[5 X P][5' X P] 


II 


+6 


-3.972(1) 





[V;. X T/p][T/p X 


I 





+16.874(1) 


— CTI 


[V,. X l/„][l/p X V^,] 


II 





-1.129(0) 





[v;, X X K] 


I 





+2.864(0) 





[v;. X K][v; X K] 


II 





-0.955(0) 





[T/^ X ^^][Vj, X Af,] 


I 





+0.201(0) 





[V^, X A^][V;, X A,,] 


II 





-0.388(0) 





[V^ X X Au] 


I 





+0.032(0) 





[V„ X ^p][V^ X aJ] 


II 





-0.096(0) 





[V^, X ^^] [V^ X Ap] 


I 





+0.032(0) 





X 4^][v;, X 


II 





-0.096(0) 





[V^ X ^,.][V^ X A^] 


I 





-0.524(0) 





[V^ X ^^][V^ X aJ] 


II 





-0.224(0) 





[T^. X 5][T^. X S] 


I 





-0.222(0) 





[T,., X 5][Tf., X S] 


II 





-0.095(0) 







I 





+0.407(0) 







II 





-1.032(0) 







I 





+1.374(0) 





X [T^t/ X 


II 





-0.458(0) 





[Tj^ii/ X r^p] \Ti_ii/ X Tiip^ 


I 





-0.127(0) 





[Tj^ll/ X r^p] \Ti_iiy X Tlyp] 


II 





+0.042(0) 





[r^tjy X r^p]['2~yii^ X r^T^] 


I 





-0.127(0) 







II 





+0.042(0) 





r,,,. X Tnn]\Tn,. X Tnv.l 


I 





—0.051(0) 





X p ] [T^ X TpT^] 


II 





+0.154(0) 





[r^tjy X Tpr^] [T^t^ X 7p,p] 


I 





-0.051(0) 





[Tf^ii/ X Tp'q'^^Ti_ii/ X ^^.p] 


II 





+0.154(0) 





[r^tjy X TpTj] [r^i^ X TpTj] 


I 





+0.276(0) 





\Tf_ii/ X jPpTj] ['Tp.t/ X T'pT7] 


II 





-0.638(0) 





[T^. X P][r^, X p] 


I 


-14/3 


-3.500(0) 





[Tm.- X P][T^.-^ X P] 


II 


-2 


-1.500(0) 






TABLE XIV. Renormalization constants of a color one trace four-fermion operator at one loop, (Ol-"")/ defined m 
Eq. (j54|). The Ttj matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The C/j"** matrix is defined in Eq. (^^. The Mij matrix represents the tadpole-improvement defined in Eq. (36). 
cti = (47r)^/12. All the greek indices are summed under the condition oi yi ^ v ^ p ^ r\. 



28 



QLatt{0) 


color trace 




/^Latt 




[A^ X V^][A^ X V^] 
[A^ X X v^,] 
[A^ X V^][A^ X K] 


I 

II 

I 







-16.069(1) 
-1.286(0) 
-2.228(0) 


+CTI 





[^M X K][Am X K] 


II 





-0.955(0) 





[Am X 1^m][Am X K] 


I 





-0.153(0) 





[Am X Vm][A^ X 
[A^ X K][A^ X V^] 
[A^ X K][A^ X V^] 


II 

I 

II 








+0.458(0) 
-0.153(0) 
+0.458(0) 







[Ap X A^][A^ X A^] 
[Am X Am] [Am X Am] 
[Ap X A„][A^ X A^] 
[Af. X A„] [A^ X A^] 










+0.110(0) 
-0.545(0) 
+0.673(0) 
-0.224(0) 


o o o o 


[P X 5'][P X 5"] 


I 





-0.077(0) 





[p X ^ilp X 5] 


II 





-0.254(0) 





[P X Tm.][P X Tm^] 


I 





-0.977(0) 





[P X Tf,,][P X T^,] 


II 





-0.419(0) 





[P X Tf,,][P X T^p] 


I 





-0.099(0) 





[P X Tf,,][P X T^p] 


II 





-0.042(0) 





[P X P][P X P] 


I 


-2 


-39.872(2) 


+2CTI 


[P X P][P X P] 


II 


+6 


+40.087(1) 






TABLE XV. Renormalization constants of a color one trace four-fermion operator at one loop, (Ol-"")/ defined m 
Eq. (g). The Vij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (^2[). 
The Cj^"** matrix is defined in Eq. (^^. The Mij matrix represents the tadpole-improvement defined in Eq. (36). 
cti = (47r)^/12. All the greek indices are summed under the condition oi ^ ^ v ^ p ^ i). 



29 





color trace 








[S X T^,][S X T^,] 
[S X T^„][S X Tp,] 
[SxT^,.][SxT,.p] 
[S X T^„][S x Tf,.p] 


I 

II 
I 

II 









-0.837(0) 
+0.279(0) 
-0.085(0) 
+0.028(0) 


o o o o 


[S X P][S X P] 


II 


+16 


-95.607(4) 


+6CT/ 


[V,i X A^][V^, X A^,] 


I 





-0.341(0) 





[V^ X A^][Vt, X A^j 


II 





+0.114(0) 





[i;, X Aj\[Vt, X A^] 


I 





-0.449(0) 





[i;, X Aj][Vt, X A^] 


II 





+0.150(0) 





[1/,, X A^,][V^, X A„] 
[v;^ X Ap] [i;, X A„] 


I 

II 






-0.096(0) 
+0.032(0) 






[Vp X AJ\[Vt, X Ap] 


I 





-0.096(0) 





[V"^ X aJ][V^ X Ap] 


II 





+0.032(0) 





[Tm-- X S][T^. X S*] 
[r„. X S][Tp, X 5] 

[r„, X P][r„, X p] 
[r„, X P][r„, X p] 


I 

II 
I 

II 





-4 
+4/3 


-0.190(0) 
+0.063(0) 
-3.000(0) 
+1.000(0) 


o o o o 


[Ap X V^] [Ap X V^] 


I 





-3.007(1) 





[A^ X V^][A^ X V^j 
[A^ X K1[A^. X K] 
[Ap X K][A^ X K] 


II 
I 

II 







+1.002(0) 
-1.909(0) 
+0.636(0) 







[Ap X Vp] [Ap X K] 
A„ X V„ A,, X V, ] 


I 

II 






+0.765(0) 
-0.255(0) 






[Ap X K] [A^ X l/p] 


I 





+0.765(0) 





[A^ X K][Aj, X 


II 





-0.255(0) 





[P X T^,][P X T^,] 


I 





-0.837(0) 





[P X T„,][P X T^,] 


II 





+0.279(0) 





[P X Tf,,][P X T^p] 
[P X Tp,][P X T^p] 


I 

II 







-0.085(0) 
+0.028(0) 






[P X P][P X P] 


II 


+16 


+111.255(2) 


-2CTI 



TABLE XVI. Renormalization constants of a color two trace four-fermion operator at one loop, {Oi)ii defined 



Eq. 



|). The Tij matrix represents an anomalous dimension defined in Eq. (Bq) and its values are given in 



Eq. 
Eq. 



69). cti 



(4^) 



The Ci," matrix is defined in Eq. (pw. The Nij matrix represents the tadpole-improvement defined in 



■/12. All the greek indices are summed under the condition oi ^ v ^ p ^ r\. 



30 



^i:,att(0) 


color trace 








[5 X S][S X 5] 
[5 X S] IS X 5] 

[SxT^,.][SxT^,] 
[S X T^„][S X T^,] 
[S X T^,.][S X T^p] 
[S X T^„][S X Tf,p] 
[S X P][5' X P] 
[S X P][5' X P] 


I 

II 
I 

II 
I 

II 
I 

II 









-6 
-6 


-0.294(0) 
-0.016(0) 
+2.000(0) 
-0.568(0) 
+0.037(0) 
-0.111(0) 
+7.343(1) 
-0.264(1) 









—cti 



[V^ X T/^][T/p X 
[V,, X X V^,] 
[v;, X K][t4 X K] 
[v;, X K][t4 X K] 

[v;, X v„][Vp X 

[14 X 14] [14 X K] 

[14 X 14] [14 X 14] 
[14 X 14] [V^, X 14] 


I 

II 
I 

II 
I 

II 
I 

II 


oooooooo 


-11.777(0) 
-0.867(0) 

-11.628(0) 
-1.318(0) 
-0.211(0) 
+0.632(0) 
-0.211(0) 
+0.632(0) 


+CT//2 



+CT1/2 







\V„ X AuWu X A,A 
[14, X 4^] [14x^1 J 
[14xA.][V4x^.] 

[I4x^.][i4x^.] 

[14 x^J[14 X A,] 

[v;, X Af,][i4 X A.] 
[v;, X A,] [14 X Ap] 
[v;, X A.] [14 X Ap] 


I 

II 
I 

II 
I 

II 
I 

II 


oooooooo 


+0.910(0) 
-0.143(0) 
+0.990(0) 
-0.385(0) 
-0.016(0) 
+0.048(0) 
-0.016(0) 
+0.048(0) 


oooooooo 



TABLE XVII. Renormalization constants of a color one trace four-fermion operator at one loop, (O^"")/ defined 
in Eq. (p4|). The Fij matrix represents an anomalous dimension defined in Eq. ( p^ and its values are given in 



Eq. 
Eq. 



The Ctj" 



matrix is defined in Eq. (|56[). The Mij matrix represents the tadpole-improvement defined in 



|66). Cti = (47r) /12. All the greek indices are summed under the condition of /i / / p 7^ r;. 



31 





color trace 




/^Latt 


Mij 


\T^. X 5][T^. X S] 


I 





+0.358(0) 





[T^u X S'lfT^. X S\ 


II 





-0.174(0) 







I 





-1.761(0) 







II 





-0.151(0) 







I 





-1.602(0) 





^Tj^i/ X r^p] [T^^ X 7/ip] 


II 





-0.627(0) 





'^[Tj.Lu X !2~'^iy][Tpi/ X Tj^ip\ 


I 





+0.042(0) 





[r^tjy X !2~'^iy][Tp,y X r^tp] 


II 





-0.125(0) 





[T^yUl^ X r^ip] [T^l/ X Tj^lJ^ 


I 





+0.042(0) 





X T/ip] [T^i/ X Tj^ij^ 


II 





-0.125(0) 







I 





+0.014(0) 





[T^^tiy X T p\^T fj^ij X T]^p^ 


II 





-0.042(0) 





[Tf^ii/ X r^p]['2~yi^< X r^T^] 


I 





+0.097(0) 







II 





-0.292(0) 





r,,,. X Tnn]\T,,,. X r^r,l 


I 





+0.026(0) 





[7pi^ X Tj_ipj \Ti_ii/ X 7p7^] 


II 





-0.077(0) 





[T^tiy X Tpr;] [T^t/ X 7/j.p] 


I 





+0.026(0) 





[Tf^ii/ X Tp'q'^^Ti_ii/ X 7/j,p] 


II 





-0.077(0) 





[T^tiy X Tp'q'^^Ti_ii/ X 7p7^] 


I 





-1.695(0) 





[7pi/ X Tprj] [T^i/ X 7pr^] 


II 





-0.349(0) 





[T;,. X P][T^, X P] 


I 


+26/3 


-6.617(1) 


+CTI 


[T;,. X P][T^, X P] 


II 


-2 


-1.915(1) 






TABLE XVIII. Renormalization constants of a color one trace four-fermion operator at one loop, (O5 ")/ defined 
in Eq. (^). The Tij matrix represents an anomalous dimension defined in Eq. ( ^^ and its values are given in 
Eq. 
Eq. 



The Cij matrix is defined in Eq. (|56|). The Mij matrix represents the tadpole-improvement defined in 
|66). cti = (47r)^/12. All the greek indices are summed under the condition oi ^ v ^ p ^ i). 



32 



^„Latt(0) 


color trace 






71 f 

Mij 


[A^ X V^][A^ X V^] 


I 





+12.329(0) 


-CTl/2 


[A^ X X V^] 


II 





-0.789(0) 





[^M X X K] 


I 





+12.480(0) 


-ctj/2 


[^M X X K] 


II 





-1.239(0) 





[^M X X K] 


I 





+0.076(0) 





[A^ X V^][A^ X 


II 





-0.229(0) 





[^M X X V^] 


I 





+0.076(0) 





[^M X K][Am X V^] 


II 





-0.229(0) 





[^M X K][Am X Vp] 


I 





-0.223(0) 





[A^ X K][Am X Vp] 


II 





+0.669(0) 





[Af, X A^] [Af, X A^] 


I 





-0.841(0) 





[Ap X aJ\ [Ap X A^j 


II 





-0.064(0) 





[Ap X A„] [Ap X A^] 


I 





-0.760(0) 





[A^ X A^][A^ X A^j 


II 





-0.306(0) 





[A^ X AJ[A^ X A^] 


I 





+0.016(0) 





[Ap X Ap] [Ap X A^] 


II 





-0.048(0) 





[Ap X A„] [A^ X aJ\ 


I 





+0.016(0) 





[Ap X A„] [Ap X A^] 


II 





-0.048(0) 





[A^ X A^][Ap, X Ap] 


I 





+0.032(0) 





[A^ X A^][Af, X Ap] 


II 





-0.096(0) 





[P X 5'][P X 5"] 


I 





-0.294(0) 





[P X ^ilp X S] 


II 





-0.016(0) 





[P X Tf,,][P X Tf,,] 


I 





+2.066(0) 





[P X Tf,,][P X T^,] 


II 





-0.765(0) 





[P X Tp,][P X T^p] 


I 





-0.014(0) 





[P X Tp,][P X T^p] 


II 





+0.042(0) 





[P X P][P X P] 


I 


-6 


+14.686(1) 


—cti 


[P X P][P X P] 


II 


-6 


-22.293(1) 






TABLE XIX. Renormalization constants of a color one trace four-fermion operator at one loop, (Ol-"")/ defined m 
Eq. (j54|). The Vij matrix represents an anomalous dimension defined in Eq. (|5^) and its values are given in Eq. (52). 
The C/j"** matrix is defined in Eq. (^^. The Mij matrix represents the tadpole-improvement defined in Eq. (36). 
Cti = (47r)^/12. All the greek indices are summed under the condition oi ^ ^ v ^ p ^ i). 



33 



^i:,att(0) 


color trace 




^Latt 




[5 X S][S X 5] 
[5 X S] IS X 5] 

[S X Tp,,][S X Tp,] 
[S X Tp^][S X Tp,\ 
[S X Tp,,][S X Tpp] 
[S X T^„][S X Tpp] 
[S X P][S X P] 
[S X P][S X F>] 


I 

II 
I 

II 
I 

II 
I 

II 









-6 

-6 


-0.286(0) 
+0.095(0) 
+0.419(0) 
-0.140(0) 
+0.042(0) 
-0.014(0) 
-4.500(0) 
+49.304(2) 










-3ct/ 


[Vp X Vp][Vp X Vp] 
[Vp X Vp][Vp X Vp] 
[Vp X V^][Vp X K] 
[Vp X V^][Vp X K] 
[Vp X Vp][Vp X 
[Vp X Vp][Vp X 
[Vp X V^][Vp X Vp] 
[Vp X V^][Vp X Vp] 
[Vp X VA[Vp X Vp] 
[Vp X V^][Vp X Vp] 


I 

II 
I 

II 
I 

II 
I 

II 
I 

II 


oooooooooo 


-2.863(0) 
+0.954(0) 
-3.413(0) 
+1.138(0) 
+0.382(0) 
-0.127(0) 
+0.382(0) 
-0.127(0) 
+0.765(0) 
-0.255(0) 


oooooooooo 


[V„ X A,j]\V„ X A„] 
[Vp X Ap][Vp X Ap] 
[Vp X A^][Vp X A,] 
[Vp X A^][Vp X A,] 
[Vp X Ap][Vp X A^] 
[V;, x^^][V;xA,] 
[Vp X A^][Vp X Ap] 
[Vp X A^][Vp X Ap] 


I 

II 
I 

II 
I 

II 
I 

II 


oooooooo 


+0.171(0) 
-0.057(0) 
+0.224(0) 
-0.075(0) 
+0.048(0) 
-0.016(0) 
+0.048(0) 
-0.016(0) 


oooooooo 



TABLE XX 
in Eq. (H). The 



Eq. 
Eq. 



Renormalization constants of a color two trace four-fermion operator at one loop: (Og"")// defined 
ij matrix represents an anomalous dimension defined in Eq. ( p^ and its values are given in 

. The Nij matrix represents the tadpole-improvement defined in 



The C,^"" matrix is defined in Eq. (|f_ 
cti = (47r)^/12. All the greek indices are summed under the condition oi ^ v ^ p ^ i). 
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color trace 








[T^. X 5][T^. X S] 


I 





+0.095(0) 





[T^u X S'lfT^. X S\ 


II 





-0.032(0) 







I 





-1.335(0) 







II 





+0.445(0) 







I 





-1.711(0) 





^Tj^i/ X r^p] [T^^ X 7/ip] 


II 





+0.570(0) 





'^[Tj.Lu X !2~'^iy][Tpi/ X Tj^ip\ 


I 





-0.201(0) 





[r^tjy X !2~'^iy][Tp,y X r^tp] 


II 





+0.067(0) 





[T^yUl^ X r^ip] [T^l/ X Tj^lJ^ 


I 





-0.201(0) 





X T/ip] [T^i/ X Tj^ij^ 


II 





+0.067(0) 





\T„„ X r„„ T„„ X r„„i 


I 





-0.445(0) 





\T'l.LU X r^p]['2~yi/y X T]^p^ 


II 





+0.148(0) 





1^1^11/ X r^p] [T^^- X r^T^] 


I 





-0.042(0) 







II 





+0.014(0) 





X rpTj]['2~yi^ X 7p7^] 


I 





-1.335(0) 





X TpT^] [T^t/ X 7p7^] 


II 





+0.445(0) 





[r„„ X P][r„„ X p] 


I 


+2 


+1.500(0) 





[r„, X P][r„, X p] 


II 


-10/3 


-15.119(1) 


+CTI 



TABLE XXI. Renormalization constants of a color two trace four-fermion operator at one loop, (Og"")// defined 
in Eq. (^). The Tij matrix represents an anomalous dimension defined in Eq. ( ^^ and its values are given in 
Eq. " " 
Eq. ( 



The Cij" matrix is defined in Eq. (|52). The Nij matrix represents the tadpole-improvement defined in 
cti = (47r)^/12. All the greek indices are summed under the condition oi ^ ^ v ^ p ^ r). 
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o 



Latt{0) 



[A^ X X V^] 

[A^ X X K] 

[^M X X K] 

[^M X X K] 

[A^ X V^][A^ X 

[^M X X Vf,] 

[A^ X X V^] 



color trace 



[Af, X X A^] 

[Af, X AJ[A^ X A,,] 
[Af, X X A^] 

[Af, X A„] [A^ X A^] 
[Af, X A^] [A^ X A„] 
[A^ X aJ\[A^, X A^j 
[A^ X A^][A^ X A^] 
fA„ X A^][A^ X A„l 



I 

II 
I 

II 
I 

II 
I 

II 



+1.504(0 
-0.501(0 
+0.955(0 
-0.318(0 
-0.382(0 
+0.127(0 
-0.382(0 
+0.127(0 



I 

II 
I 

II 
I 

II 
I 

II 



-0.673(0 
+0.224(0 
-1.014(0 
+0.338(0 
-0.143(0 
+0.048(0 
-0.143(0 
+0.048(0 



[P X S][P X S] 
[P X S] [P X S] 



[P X T^,][P X T^p] 
[P X Tf,,][P X T^p] 
[P X P][P X P] 
[P X P][P X P] 



I 

II 
I 

II 
I 

II 

I 

II 










-6 

-6 



-0.286(0 
+0.095(0 
+0.419(0 
-0.140(0 
+0.042(0 
-0.014(0 
-4.500(0 
-54.128(1 











+CTI 



TABLE XXII. Renormalization constants of a color two trace four-fermion operator at one loop, (Of""")// defined 
in Eq. (p4|). The Fij matrix represents an anomalous dimension defined in Eq. ( p^ and its values are given in 
Eq. " 
Eq. 



The Cij" matrix is defined in Eq. (jS^)- The Nij matrix represents the tadpole- improvement defined in 
|69[). cti = (47r)^/12. All the greek indices are summed under the condition of fi ^ u ^ p ^ rj. 
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